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In this work we present the preliminary results of a complete equation of state (EoS) for core¬ 
collapse supernova simulations. We treat uniform matter made of nucleons using the the quark- 
meson coupling (QMC) model. We show a table with a variety of thermodynamic quantities, which 
covers the proton fraction range Yp = 0—0.65 with the linear grid spacing AY), = 0.01 (66 points) and 
the density range ps ~ 10^“^ — 10^®g.cm“® with the logarithmic grid spacing Alogio(pfl/[g.cm“^]) = 
0.1 (21 points). This preliminary study is performed at zero temperature and our results are 
compared with the widely used EoS already available in the literature. 


I. INTRODUCTION 

Although the theory related to supernova (SN) has 
made a remarkable progress in the past decade, there are 
many questions still to be answered. The catastrophic 
infall of the core of a massive star, reversed to trigger the 
powerful ejection of the stellar mantle and envelope in 
a supernova explosion, was identified as the crucial role 
in the synthesis of heavy elements. But when, why and 
how it happens is a fundamental problem of stellar as¬ 
trophysics that remains to be explained. The implosion 
of stellar cores was als^roposed as part of the scenario 
of the stellar death [I|, |J. 

There are no doubts that supernovae explosions are an 
unique phenomenon in nature and an excellent labora¬ 
tory to test extreme physics conditions. Unfortunately 
SN tends to happen once in a century per galaxy which 
makes the study of these great labs very difficult. Due to 
this observational difficulty, simulations of core-collapse 
supernova have played an important role in the study of 
supernovae explosions and their possible remnants. The 
equation of state of nuclear physics is a fundamental in¬ 
gredient in the simulation of SN explosion. Simulations of 
core-collapse supernovae have to be fed with a wide range 
of thermodynamic conditions. Moreover, extremely high 
density and temperature may be achieved when black 
holes are formed by failed supernovae. To date, the tem¬ 
perature is believed to vary from zero to more than 100 
MeV, densities from 10® to more than 10^® g.cm“® and 
the proton fraction up to about 0.6. Fulfilling these con¬ 
ditions makes the construction of a complete EoS a very 
hard work, mainly at low densities, where a variety of 
sub-structures and light clusterization are possible. For 
these reasons, there are only few complete EoS available 
in the literature. The most commonlyused EoS are those 
of Lattimer and Swesty Q , Shen [1, and Hempel @ . 

The Lattimer EoS is based on a compressible liquid 
drop model with a Skyrme force for nucleon interactions. 
The Shen EoS from 1998 was the first equation of 
state for supernova simulations using a relativistic nu¬ 
clear model. The upgrade to Shen’s work was published 
in 2011 with more points in the table and the inclu¬ 
sion of the A hyperons. Both works developed by Shen 
were constructed with the relativistic mean field Walecka 
model 0 using the TMl Q parameterization. Hempel’s 


EoS [1] is based on the TMl, TMA 0 and FSUgold 
0 parameterizations of the Walecka model and use 
the nuclear statistical equilibrium model of Hempel and 
Schaffner-Bielich [Tlj, which includes excluded volume 
effects. 

An interesting analysis of the different parameteriza¬ 
tions of the relativistic-mean-field (RMF) models was 
made by Dutra et al [l^, where the authors analyzed 
263 different RMF models under several constrains re¬ 
lated to the symmetric nuclear matter (SNM) and pure 
neutron matter (PNM). The TMl parameterization used 
by Shen failed under six of this constraints, the TMA 
used by Hempel failed in five, and the FSUgold, also 
used by Hempel failed in one constraint. More details of 
this analysis in (l^ . 

The works of Hempel, Lattimer and Shen were suc¬ 
cessful and very useful in many calculations in the last 
decade [IM3 , but there are some simulations of SN in 
which the supernova does not explode 0 . It is believed 
that these failures are due to some problems with the 
nuclear EoS. 

In this work we present our preliminary results for 
the construction of an EoS grid for core-collapse super¬ 
nova simulations, with the quark-meson coupling (QMC) 
model 0. 

In the QMC model, nuclear matter is described as a 
system of non-overlapping MIT bags [2l| that interact 
through the exchange of scalar and vector meson fields. 
Many applications and extensions of the model have been 
made in the last years [2^ - 10 . It is found that the EoS for 
infinite nuclear matter at zero temperature derived from 
the QMC model is softer than the one obtained with the 
Walecka model. This might be a problem if one wants to 
describe very massive compact objects [13,0, but as far 
as SN simulations are concerned, it is worth testing it be¬ 
cause apart from numerical accordance, one can interpret 
that starting from quark degrees of freedom is an advan¬ 
tage on the physical meaning underneath. Moreover, the 
effective nucleon mass obtained with the QMC model lies 
in the range 0.7-0.8 of the free nucleon mass, which agrees 
with the results derived from the non-relativistic analysis 
of scattering of neutrons from lead nuclei 0 | and is larger 
in comparison with the effective mass obtained with some 
of the different parameterizations of the Walecka models. 
A low effective mass at saturation can be a problem when 
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Model 

M 

(MeV) 

rriq 

(MeV) 

rria 

(MeV) 

niuj 

(MeV) 

rrip 

(MeV) 

9p 

g^j 

9<y 

(MeV) 

NLT 

DDF 

QMC 

939.0 

5.5 

550 

783 

770 

8.6510 

8.9817 

5.9810* 

210.85 

no 

no 

NL3 

939.0 

- 

508.194 

782.501 

763 

8.9480 

12.868 

10.217 

- 

yes 

no 

GMl 

939.0 

- 

550 

783 

770 

8.1945 

10.608 

9.5684 

- 

yes 

no 

TMl 

938.0 

- 

511.197 

783 

770 

4.6321 

12.613 

10.028 

- 

yes 

no 

FSUgold 

939.0 

- 

491.5 

782.5 

763 

11.7673 

14.301 

10.592 

- 

yes 

no 

TMAt 

938.9 

- 

519.151 

781.95 

768.1 

3.800 

12.842 

10.055 

- 

yes 

no 

TW 

939.0 

- 

550 

783 

763 

7.32196“ 

13.2901“ 

10.728“ 

- 

no 

yes 


TABLE I: Parameters used in the QMC model and different Walecka model parameterizations. In the first line we present the 
parameters used in the present work with the QMC model. The NL3, GMl and TW parameterizations of the Walecka model 
are used here for a comparison of bulk matter properties. TMl is the parameterization used in Shen’s work. The FSUgold and 
TMA were used in Hempel’s work. *g‘^ is the quark-meson coupling in the QMC model. ** values taken at saturation for the 
TW parameterization. TMA^: the coupling parameters gt of the set TMA are chosen to be mass-number dependent such that 
gi = ai + with at and bi being constants [^; for infinite matter as in the stellar matter, one has an infinite nucleus, 

and then the limit A i—>■ infinity is taken so that gi = ai. NLT= Nonlinear terms. DDP= Density dependent parameters.) 


hyperons are included in the calculation, as discussed in 
the next section. 

All the few works for the purpose of the SN simula¬ 
tions that are based on relativistic models, use different 
parameterizations of the Walecka model, in which the nu¬ 
cleons interact between each other through the exchange 
of mesons. We believe that with the quarks degree of 
freedom present in the QMC model, a more fundamen¬ 
tal physics lacking in the already used models, can be 
tested and possibly contribute for the SN simulations to 
explode. 

Another possible use of this preliminary EoS obtained 
with the QMC model is the study of the cooling of com¬ 
pact stars, which serve as an important window on the 
properties of super-dense matter and neutron star struc- 
ture, and is very sensitive to the nuclear equation of state 

This paper is structured as follows: in the second sec¬ 
tion we present a review of the quark-meson coupling 
model. In the third section we present our results for the 
EoS. The last section is reserved for the final remarks 
and conclusions. 


II. THE QUARK-MESON COUPLING MODEL 


malized ground state for a quark in the bag is given by 
exp(-ieg^t/i?Ar) 


Xq 


where 


and. 
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( 4 ) 


+RNm* 

with the normalization factor given by 

K" = ‘^RNJoixq) - 1) + RNm*q/2] /xl , (5) 


where 0,^ = -b {Rn m*)^, m*q =m°- g% a, Rn is 

the bag radius of nucleon N and Xq is the quark spinor. 
The bag eigenvalue for nucleon N, Xq^^ , is determined by 
the boundary condition at the bag surface 


In the QMC model, the nucleon in nuclear medium is 
assumed to be a static spherical MIT bag in which quarks 
interact with the scalar (a) and vector (oj, p) fields, and 
those are treated as classical fields in the mean field ap¬ 
proximation (MFA) [ 23 . The quark field, ipqN^ inside the 
bag then satisfies the equation of motion: 


+ ^gpXzP03l° fPqNix)=0 


q = u,d (1) 


where rriq is the current quark mass, and g'^, g^ and 
denote the quark-meson coupling constants. The nor- 


J 0 jv(^ 9 n) l^qN JxwiXqN) ■ 


( 6 ) 


The energy of a static bag describing nucleon N con¬ 
sisting of three quarks in ground state is expressed as 




n. 


QN 


Rn 


Zn 

Rn 


- TT Bm , (7) 


where Zjy is a parameter which accounts for zero-point 
motion of nucleon N and is the bag constant. The 
set of parameters used in the present work is determined 
by enforcing stability of the nucleon (here, the “bag”), 
much like in [33| , so there is a single value for proton and 
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neutron masses. The effective mass of a nucleon bag at 
rest is taken to be 

The equilibrium condition for the bag is obtained by 
minimizing the effective mass, with respect to the 
bag radius 


dM% 

dR*^ 


N = p,n, 


( 8 ) 


By fixing the bag radius Rn = 0.6 fm and the bare nu¬ 
cleon mass M = 939 MeV the unknowns = 4.0050668 
and = 210.85MeV are then obtained. Furthermore, 
the desired values of B/A = tjp — M = —15.7 MeV at 
saturation n = no = 0.15 fm“^, are achieved by setting 
= 5.9810, = 8.9817, where g^ = igl and gp = 3g«. 

All the parameters used in this work are shown in Table 

III 

The total energy density of the nuclear matter reads 

^ + irUpPoa 

-1 pkN 

+ E— / k^dk[k^-M*^^]^/^, (9) 

N ^ do 


and the pressure is. 



k^dk/[k^ - M^2]1/2^ 


( 10 ) 


The vector mean field ojq and po 3 are determined 
through 


gujij^p A nri^ riyj) 

Wo =-5-, Po3 =-5-: 

—T ^ nry-, ^ 


mt 




(11) 


where 


ns = 


2fc3 


'TV 


37r^ 

N 


, N = p,n. 


( 12 ) 


is the baryon density. 

Finally, the mean field cr is fixed by imposing that 



(13) 


It is always important to check the behavior of the 
models in the symmetric nuclear matter at saturation 
density and zero temperature, i.e., the bulk nuclear mat¬ 
ter properties. A comprehensive work in this direction is 
[T^ . but the QMC model was not analyzed. Therefore, 
we compare the QMC model with two parameterization 
of the well known Walecka-type models, namely: GMl 
[S^ and N L3 and the density dependent parameter 
model TW [^. In this work we have chosen GMI for be¬ 
ing a parameterization which gives a good value for the 
effective mass, NL3 because it is a very common stan¬ 
dard parameterization and TW because it is a very good 



FIG. 1: The effective mass of the QMC model and three 
Walecka parameterizations. 



FIG. 2: Equation of state of the QMC model and three 
Walecka model parameterizations. 


Model 

B/A 

no 

M*/M 

^sym 

K 


(MeV) (fm-3) 


(MeV) (MeV) 

QMC 

-15.7 

0.150 

0.77 

34.5 

295 

NL3 

-16.2 

0.148 

0.60 

37.4 

272 

GMl 

-16.3 

0.153 

0.70 

32.5 

300 

TMl 

-16.3 

0.145 

0.63 

36.8 

281 

FSUgold 

-16.3 

0.148 

0.62 

32.6 

230 

TMA 

-16.0 

0.147 

0.63 

30.7 

318 

TW 

-16.2 

0.153 

0.56 

32.6 

240 


TABLE II: Nuclear matter bulk properties obtained with 
the QMC model, two different parameterizations of Walecka 
model and one density dependent model we use in this paper, 
and the four parameterizations used in the works of Hempel 
and Shen. All quantities are taken at saturation. 




















4 



FIG. 3: Binding energy of the nucleons as function of the 
baryon density with the proton fractions Yp =0, 0.1, 0.3, 0.5 
and 0.6. 



FIG. 5: The relation between pressure and pB for the 
Yp =0.01, 0.1, 0.3, 0.5 and 0.6 protons fractions. 
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FIG. 4: Compression modulus as function of the baryon den¬ 
sity with the proton fractions Yp =0.1, 0.3, 0.5 and 0.6. 



density dependent parameterization according to [l^ . In 
[l^ . it was found that GMl failed under six constrains 
related to the symmetric nuclear matter (SNM) and pure 
neutron matter (PNM), NL3 failed under nine and TW 
did not satisfy only one constraint. 

In Table U we can see the different parameters here 
analyzed and the ones used in the works of Shen an 
and Hempel Q. 

In Table im we show the properties of nuclear matter at 
saturation density and zero temperature. The first col¬ 
umn shows the binding energy per nucleon, the second 
shows the baryonic density of the nucleons at the satura¬ 
tion point, the third the effective mass relative to the free 
nucleon mass, the symmetry energy and the compression 
modulus of the nucleons. 

In Fig. [T] we show the relation of the effective mass as 
function of the baryonic density where we can see that 


the effective mass for the QMC model has always a larger 
value than the other models, which is an important char¬ 
acteristic if hyperons are to be included. As can be seen 
in [s^ the effective mass of some of the Walecka-type 
models tend to zero at still low densities when all the 
baryons of the octet are included. For this reason, param- 
eterizations with larger effective masses were proposed by 
Glendenning [s^ (GMl, GL, etc) with the specific pur¬ 
pose of applications to stellar matter. Note that in both, 
figure and table, we can see that M*/M = 0.77 at the 
saturation density, which agrees with the results derived 
from the non-relativistic analysis of scattering of neu¬ 
trons from lead nuclei [2^, and that this result is larger 
in comparison with the effective mass from many of the 
Walecka-type models. 

In Fig. [2] we show the pressure versus energy density 
relation for infinite nuclear matter at T = 0, where we 
see that the curve for QMC is similar to the one for TW, 
and softer than the GMl and NL3. As mentioned in 
the Introduction, this can be a deficiency for neutron 
star studies, but those effects on supernovae are still not 
analyzed. 


III. RESULTS AND DISCUSSIONS 

In this work we construct an EoS table covering a wide 
range of proton fraction Yp and baryon density ub- We 
show in the paper some results concerning the properties 
of matter with the QMC model and indicate the website 
from where the full table can be download. 

We show in Fig. |3]the binding energy of homogeneous 
nuclear matter at zero temperature as a function of the 
baryon density for different proton fractions. For pure 
neutron matter and low proton fractions, there are no 
binding states, as expected. This result is in good agree¬ 
ment with Walecka Q and Shen [^. In Fig. |4]we show 
the compression modulus versus the baryonic density. 
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FIG. 6: Some fixed values of proton fraction for chemical 
potentials as function of the baryonic density. The continnous 
line represents fj.p and the dashed line represents /in. 


where the expression for K is 


K = 9 rig 




drig 



(14) 


ns is the baryon number density and e the total energy 
density of the nucleons. Here we see that near the satura¬ 
tion density the compression modulus has the same value 
independent of the proton fraction and, as the density in¬ 
creases, the compression modulus is different for each Yp. 

In Fig. Othe pressure p as a function of ps is displayed. 
The baryon number density is related to the baryon mass 
density as pB = where = 931.49432 MeV is 

the atomic mass unit. We note that the pressure varies 
more with the Yp for low values of p at lower densities, 
this result is also in good agreement with Shen’s work 

i. 

In Fig. [S] we show the proton and neutron chemical 
potentials, Pp and /j,„, as function of the baryon density 
for the proton fractions Yp = 0.1, 0.3, 0.5 and 0.6. In 
these curves we see that for = 0.1 the chemical poten¬ 
tial of the neutron is bigger than the one of the proton. 
As the proton fraction gets bigger, the curves approach 
each other, until they are the same in = 0.5, and for 
Yp = 0.6 the chemical potential of the proton is bigger 
than the one of the neutron. This is an obvious result, but 
as the chemical potentials are very important quantities 
in the EoS tables, they are also presented graphically. 

Once bulk nuclear matter properties are shown to be¬ 
have as expected and present some important differences 
as compared with the other works, we proceed toward 
building a preliminary EoS table with the QMC model, 
for homogeneous matter and zero temperature, which is 
available on the Web at 


http://debora.fsc.ufsc.br/eos_qmc.tO, 

In Table Hill we show the thermodynamic quantities 
described as in Si 

1. Temperature: T[MeV]. 

2. Logarithm of baryon mass density: 

logio(PB)[g-cm-3]. 

3. Baryon number density: ub [fm“^]. 


4. Proton fraction: Yp. 

The proton fraction Yp of uniform matter made of 
protons and neutrons, is defined by 


-I- rip 

where tip and rin are the number density of protons 
and neutrons, respectively. 

5. Free energy per baryon: F [MeV]. 

The free energy per baryon reads, 

f = e-Ts. 




















































6 


T 

(MeV) 

logio (ps) 
(g.cm-®) 

UB 

(fm-®) 

Yp 

F Ei-at S 

(MeV) (MeV) (kp 

M*jv 
) (MeV) 


Xp 

p 

(MeV fm-®) 

l^p 

(MeV) (MeV) 

0 

14.0 

0.0602 

0 

4.890 

12.40 

0 

838.9 

1 

0 

0.2371 

23.43 

-68.43 

0 

14.1 

0.0758 

0 

6.090 

13.60 

0 

816.9 

1 

0 

0.5103 

31.20 

-82.21 

0 

14.2 

0.0954 

0 

8.133 

15.64 

0 

791.2 

1 

0 

1.0890 

42.68 

-97.63 

0 

14.3 

0.1201 

0 

11.56 

19.07 

0 

761.5 

1 

0 

2.2780 

59.66 

-114.3 

0 

14.4 

0.1512 

0 

17.19 

24.70 

0 

728.8 

1 

0 

4.6550 

84.64 

-131.4 


TABLE III: EoS table at T = 0. It covers the proton fraction range Yp — 0 — 0.65 with the linear grid spacing AYp = 0.01 (66 
points), and the density range ps ~ — 10*^® g cm“® with the logarithmic grid spacing Alogio (ps /[g cm“®]) = 0.1 (21 

points). This table is available in the website . An exerpt of it is shown here for guidance. 


This work is for zero temperature only, hence / = e. 
The free energy per baryon is defined relative to the 
nucleon mass as, 

F = — - M = B/A. 

riB 

6. Internal energy per baryon:i?i„( [MeV]. 

The internal energy per baryon is defined relative 
to the atomic mass unit m„ = 931.49432 MeV as 

e 

^int — 

riB 

7. Entropy per baryon: S'[/cb]. 

The case we work here is for zero temperature, 
therefore S = 0. 

8. Effective nucleon mass: [MeV]. 

The effective nucleon mass is obtained in the QMC 
model for uniform matter with the relation = 
where N = p, n, and the bag energy is ob¬ 
tained through the Eq. 0 

9. Free neutron fraction: 

10. Free proton fraction: Xp. 

11. Pressure: p [MeV.fm“^]. 

The pressure is calculated from Eq. [TUI 

12. Chemical potential of the neutron: /i„ [MeV]. 

For the zero temperature case, the chemical poten¬ 
tial of the neutron relative to the free nucleon mass 
M is calculate through the following equation 

^J■n = [kn + -|- PcjOJo — ^P03 “ M. 
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IV. CONCLUSION AND FUTURE WORKS 

In this work we have used the QMC model for the first 
time for the construction of a preliminary EoS that in 
the future can be useful for the studies involving cooling 
of neutron stars and supernova simulations. We believe 
that with the quarks degree of freedom present in the 
QMC model, the EoS can contribute with part of the 
physics that lack for SN simulations to explode. 

The next step of the present work, already under de¬ 
velopment, is the computation of the EoS grid at finite 
temperature, which is essential for the supernova simu¬ 
lations. 

Thereafter we will study the very low density regions, 
where nuclear matter is no longer uniform. This will be 
done with the pasta phase approach [s^, |4^ . We believe 
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